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INTRODUCTION 

The  failure  diameter  of  an  explosive  composition  is  an  important 
characteristic  which,  to  some  extent,  determines  the  composition's 
usefulness  from  an  application's  viewpoint.  The  failure  diameter  is 
that  dimension  through  which  a steady  detonation  cannot  be  propagated. 
One  class  of  explosives  which  has  evoked  a great  deal  of  interest  in 
recent  years  has  been  labeled  "non-ideal"  explosives.  This  class 
of  explosives  is  not  thought  to  react  completely  within  the  steady 
state  detonation  wave.  For  this  discussion  non-ideal  explosives  are 
compositions  containing  an  ideal  explosive  (RDX,  HMX) , aluminum,  and 
an  oxidizer  (ammonium  nitrate,  ammonium  perchlorate) . The  failure 
diameters  of  non-ideal  explosives  are  usually  much  greater  than  those 
of  ideal  explosives.  The  usual  method  of  determining  failure  diameter 
is  to  fire  a series  of  cylindrical  test  charges  of  different  diameter. 
This  is  an  expensive  and  time  consuming  process,  especially  if  one  is 
developing  a new  composition.  The  ultimate  goal  of  this  work  is  to 
predict  the  failure  diameter  of  an  explosive  composition  on  the  basis 
of  the  physical  properties  of  its  constituents. 

To  be  able  to  predict  failure  diameter  in  this  manner,  the 
physical  properties  related  to  failure  diameter  must  be  determined. 

Two  properties  which  seem  clearly  related  to  detonation  failure 
are  the  speed  of  sound,  c,  in  the  steady  detonation  wave  and  the 
reaction  zone  length,  1.  Rarefaction  waves  entering  from  the 
cylindrical  edges  of  the  steady  detonation  front  can  certainly  quench 
or  slow  reaction  if  they  penetrate  deeply  enough  into  the  steady 
wave.  They  can  also  cause  curvature  of  the  detonation  front.  The 
depth  of  penetration  will  depend  on  the  sound  speed  and  the  time 
available  for  propagation  into  the  steady  wave  regime.  The  latter 
time  will  in  turn  depend  on  reaction  zone  length.  It  is  rather 
interesting  that  in  adding  aluminum  to  explosives  both  of  these 
quantities  may  be  expected  to  increase.  The  sound  speed  which 
ordinarily  decreases  when  small  quantities  of  aluminum  are  added  to 
explosives  will  probably  increase  (due  to  the  same  aluminum  presence) 
in  the  dense  regions  behind  a detonation  front.  In  addition,  the 
slower  energy  release  rate  associated  with  reacting  aluminum  should 
extend  the  length  of  the  steady  detonation  wave  (reaction  zone) . Such 
explosives  should,  therefore,  be  expected  to  have  larger  failure 
diameters  than  their  non-aluminized  counterparts. 

diameter  effect 

The  phenomenon  known  as  diameter  effect  in  condensed  explosives 
is  no  more  than  the  experimental  fact  that  the  detonation  velocity,  D, 
of  a cylindrical  charge  decreases  as  the  diameter  of  the  charge 
decreases.  This  relationship  between  detonation  velocity  and  diameter 
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has  been  studied  theoretically  by  Jones, ^ Eyring,  Powell,  Duffey, 
and  Parlin,2  wood  and  Kirkwood,^  Sichel,^  and  Dubnov,^  to  mention  a 
few.  (Actually,  Wood  and  Kirkwood  were  concerned  with  radius  of 
curvature  of  the  detonation  wave  and  not  the  diameter) , The 
relationships  which  have  been  derived  by  these  researchers  can  be  put 
into  the  form. 


where  is  the  detonation  velocity  at  infinite  diameter,  I is 
the  reaction  zone  length,  R is  the  radius  of  the  charge,  and  A is 
a constant  (which  is  dependent  on  the  model  used  to  derive  it) . This 
equation  could  be  used  to  predict  failure  diameter  (1)  if  £ was  a 
constant  for  varying  R or  (2)  if  a relationship  between  £ and  R 
could  be  obtained.  If  condition  (1)  were  true,  the  problem  would  be 
trivial.  However,  £ is  not  constant  as  will  be  shown  below. 

Condition  (2)  is  a possibility  because  £ does  vary  with  R. 

MODIFICATION  OF  D VERSUS  1/R  EQUATION 

The  reaction  zone  length  of  an  explosive  is  related  to  the 
sound  speed  in  the  reaction  zone  and  the  reaction  time,  t , by 

T 

£ = / c(t)  dt.  (2) 

o 


At  the  present  time  c(t)  is  an  unknown  function  which  depends  on 
the  reaction  kinetics  of  the  explosive.  In  a steady  detonation,  c(t) 
is  thought  to  be  a monotonically  decreasing  function  of  time,  thus 
Equation  (2)  can  always  be  rewritten  in  the  form 

£ = c T (3) 

T 

where  c^  = ^ / c(t)  dt,  and  should  satisfy  ^ c ^ ^CJ'  ^f 

o 

sound  speed  at  the  detonation  front,  and  c^j  is  the  sound  speed  at 
the  end  of  the  reaction  zone.  Substituting  Eq. (3)  into  Eq.  (1)  yields 


1.  H.  Jones,  Proc.  Roy.  Soc.  (London) , A189,  415  (1947). 

2.  Eyring,  Powell,  Duffey,  and  Parlin,  Chem.  Revs.,  45 , 69,  (1949). 

3.  Wood  and  Kirkwood,  Jour.  Chem.  Phys.")  22,  1920-24  (1954). 

4.  M.  Sichel,  AlAA  Jour.,  4,  2^4-72  (1966). 

5.  L.  V.  Dubnov,  Russ.  Jour.  Phys.  Chem.,  34  (10),  1124-25  (1960). 
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To  proceed  further  requires  a set  of  data  on  an  explosive  which 
gives  D and  x as  a function  of  R for  various  R's  including  the 
failure  radius,  Rf.  Dremin,  Saviov,  Trofimov,  and  Shvedov®  have 
published  such  data  for  cast  TNT  at  densities  of  1.62  and  1.60  g/cc . 
Table  1 lists  this  data  (Table  2 of  reference  6) . The  difference 
in  density  (and  failure  radius)  is  due  to  different  methods  of 

casting.  " The  trotyl  [TNT]  castings  of  the  first  type  were 

obtained  by  cooling  the  melt  with  the  continuous  stirring.  They  had 
a homogeneous  fine-crystalline  structure  throughout  the  entire  cross 
section  of  the  charge.  The  density  of  such  castings,  determined  by 
analytic  suspension  in  water,  amounted  to  1.62  g/cm^ . Castings  of 
the  second  type  were  manufactured  similarly,  but  the  melt  was  cooled 
without  stirring.  As  a result  castings  were  obtained  with  large 
crystals,  extending  from  the  periphery  to  the  center,  and  their 
density  was  1.6  g/cm^  on  the  average...  PCJ  and  ucj  in  Table  1 

are  the  pressure  and  particle  velocity  at  the  end  of  the  reaction 
zone.  The  detonation  parameters  in  Table  1 were  obtained  by  Dremin 
as  follows.  D was  determined  by  streak  camera  measurements,  ucj  and 
T were  obtained  from  the  break  point  in  the  u(t)  profile  obtained 
from  electromagnetic  velocity  gages  located  on  the  axis  of  the  charge. 
Pcj  was  obtained  from  Pcj  = Po  ^ was  obtained  from  l=i  (D-1.2  5 

Ucj)-*  Since  there  is  only  a 1%  difference  in  density,  one  value  of 
Dq  (7.04  mm/ysec)  will  be  used  for  these  two  densities. 

It  is  obvious  from  Table  1 that  i is  not  a constant  for  varying 
R.  The  next  question  is:  Is  A of  Equation  (1)  a constant?  Solving 
Equation  (1)  for  A one  obtains 

A=(l-|^)|-  (5) 

For  1.62  g/cc  TNT,  A varies  from  0.071  at  R = 30  mm  to  0.461  at 
R = 8 mm. 

Up  to  this  point,  only  the  variation  of  ^ with  1 and  R has 

been  considered.  The  data  of  Table  1 can  also  be  used  to  look 
at  the  variation  of  Pjj  with  I and  R.  An  equation  similar  to 
Equation  (1)  or  Equation  (4)  involving  Pgj  can  be  obtained  using 
the  approach  of  Dubnov.^  The  Infinite  Diameter  C-J  pressure, 

PCJo»  of  =in  explosive  is  defined^  by 

Pcjo  = G(Po)  Q (Y  - 1)  (6) 

♦Russian  researches  have  noted  that  the  extrapolated  particle  velocity 
at  the  detonation  front  using  the  EMV  gage  is  usually  1.5-1. 6 u^j. 

The  particle  velocity-time  profiles  obtained  are  essentially  linear. 

Thus  t.hey  use  the  average  value  of  the  particle  velocity  in  the  reaction 
zone  (1.25  u^j)  to  compute  1.  This  approach,  however,  ignores  any 
possibility  of  an  induction  zone  behind  the  detonation  front. 


5.  L.  V.  Dubnov,  Russ.  Jour.  Phys.  Chem.,  34 , (10),  1124-25  (1960). 

6.  Dremin,  Saviov,  Trofimov,  and  Shvedov,  Detonation  Waves  in 
Condensed  Media,  (English  Translation) . 
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Table  1 

Detonation  Parameters  of  Cast  TNT  versus  Charge  Radius* 


R 

D 

'^CJ 

^CJ 

T 

i 

mm 

mm/y sec 

mm/y sec 

Kbar 

y sec 

mm 

Pq  = 1.62 

g/cc 

30 

6.9^ 

1.62 

183.2 

0.26 

1.27 

20 

6.95 

1.60 

180.1 

0.28 

1.36 

14 

6.89 

1.52 

169.5 

0.30 

1.47 

12.5 

6.84 

1.36 

150.7 

0.31 

1.57 

11 

6.75 

1.31 

143.2 

0.30 

1.51 

10 

6.64 

1.24 

133.3 

0.34 

1.71 

9 

6.52 

1.20 

126.7 

0.33 

1.64 

8** 

6.36 

1.10 

113.3 

0.36 

1.77 

o 

II 

o 

Q. 

g/cc 

30 

6.85 

1.60 

175.5 

0.30 

1.43 

20 

6.74 

• 1.49 

159.6 

0.315 

1.51 

17.5 

6.70 

1.43 

152.3 

0.31 

1.50 

15 

6.37 

1.39 

140.8 

0.34 

1.55 

13.75**6.20 

1.20 

119.0 

0.35 

1.67 

* From  Table  2 of  Reference  6 

**  Failure  radius  of  the  charge 
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where  G(pq)  is  a function  of  the  initial  density  of  the  explosive,  Q 
is  the  total  heat  of  the  explosion,  and  y is  the  index  of  polytropy 
of  the  detonation  products.*  Taking  into  account  the  losses  due  to 
the  lateral  spreading  of  the  products  away  from  the  reaction  zone. 
Equation  (6)  is  rewritten  as 

Pqj  = G(Pq)  Q (y  -l)n  (7) 

where  P^j  is  the  actual  pressure  and  n is  a factor  characterising 
that  part  of  Q that  goes  into  the  shock  wave.  Thus 


Taking  the  square  root  of  Equation  (8)  one  obtains 


V 

j = . 


Equation  (9)  is  similar  to  Equation  (4)  of  reference  (5) . From 
reference  (5)  , ri  is  the  ratio  of  unaffected  volume  to  the  total 
volume  in  the  reaction  zone, and 

, . vl/2 


(-fe) 


^ ■ 2R 


Since  the  coefficient  is  dependent  on  the  approach.  Equation  (10)  will 
be  generalized  to 

1/2 


('^'1 
\ ^CJo  ) 


= 1 - 


Pcjo  (in  k bar)  is  obtained  from 

^CJo  “ Po  '^CJo  °o  (i2) 

where  Dq  = 7.04  mm/ysec  for  both  densities.**  u<-jq  for  the  TNT 
in  Table  1 is  obtained  by  using  y = 3 in 

^CJo  " Y + 1 ' 

*The  assumption  is  made  here  that  detonation  products  are  describable 
by  an  isentropic  equation  of  the  form  Pv^  = constant. 

**This  value  was  obtained  by  fitting  D vs  1/R. 

5.  L.  V.  Dubnov,  Russ.  Jour.  Phys.  Chem.,  34  (10),  1124-25  (1960). 
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with  the  result  that  u^jq  = 1.7  6 miri/usec.  Thus  Pqjo  ~ 2 00.7  Kbar  for 
Pq  = 1.62  and  Pcjo  ~ 198.2  Kbar  for  Pq  = 1.60  g/cc.  Table  2 lists 
A,  D/Dq,  a',  and  (Pcj/^CJo^ 1.62  g/cc  TNT,  A varies  by  a 
factor  of  2.39  while  A'  varies  by  a factor  of  1.46;  for  1.60  g/cc, 

A varies  by  1.73  while  A'  varies  by  1.50. 


Changing  from  D to  P^j  does  not  seem  to  matter  as  fer  as  the  ; 

terms  A and  A'  are  concerned.  However,  close  examination  of  the  data 
in  Table  2 shows  that  Pqj  is  more  sensitive  to  R than  D is.  For  Pq  ; 

= 1.62  g/cc  TNT  Pqj  at  failure  diameter  is  reduced  to  56%  of  its  |: 

infinite  diameter  value  while  D is  reduced  to  only  90%;  for  Pq  = 1.80  , 

g/cc  TNT  PcJ  1®  reduced  to  60%  while  D is  reduced  to  88%.  ' 

_ In  Equation  (4),  c can  be  related  to  Pqj  as  follows.  The  quantity  :{ 

c lies  between  Cf  and  Cpj  and,  it  is  reasonable  to  assume  that  c can  jij 

be  written  as  a linear  function  of  Cf  and  Cqj.  j 


The  quantity  cqj  is  given  by 

°CJ  ^ ° “ ^CJ  ° 

Since  it  follows  from  Equations  (12)  and  (13)  that  D is  proportional  to 
/pQj,  then  the  result  c^j  ® also  follows  from  Equation  (14) . 

The  quantity  Cf  can  be  closely  approximated  by  terms  linear  in  D and 
uf,  i.e.  Cf  = D +buf  where  b is  a constant  approximately  equal  to 
two.  Since  D=  /P^j,  and,  since  it  has  been  observed  that  uf  =1.5  ucjr 
it  follows  th£  Cf  is  also  approximately  proportional  to  /Pqj . 
Accordingly,  c is  written  in  the  form 


c = A" 


/p  /p 

CJo 


where  Equation  (15)  has  been  normalized  to  Pcjo* 

The  next  question  is:  Can  T be  related  to  Ppj* 

Kasley'  have  proposed  a critical  energy  fluence,  E^, 
initiation  of  explosives 

Po 


Wal)cer  and 
for  shoc)c 


where  P is  the  pressure,  t is  the  duration  of  constant  pressure  shock 
pulse,  and  Ug  is  the  shock  velocity  in  the. explosive  at  pressure  P. 
Does  an  equation  similar  to  Equation  (12)  hold  for  the  steady 
detonation  regime?  Three  possible  equations  were  investigated; 

^ 2 _ 


Pf  2 


7.  Walker  and  Wasley,  Comb  and  Flame,  22,  53-58  (1974^. 
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= 1.62  g/cc 


30 

0.991 

0.213 

0.955 

1.063 

20 

0.987 

0.191 

0.947 

0.779 

14 

0.979 

0.200 

0.919 

0.  771 

12.5 

0.972 

0.223 

0.866 

1.067 

11 

0.959 

0.299 

0.845 

: . 129 

10 

0.943 

0.333 

0.815 

1.082 

9 

0.926 

0.406 

0.795 

1.125 

8 

0.903 

0.438 

0.751 

1.125 

= 1.60 

g/cc 

30 

0.973 

0.566 

0.941 

1.238 

20 

0.957 

0.569 

0.897 

1.364 

17.5 

0.951 

0.571 

0.877 

1.435 

15 

0.905 

0.919 

0.842 

1.529 

13.75 

0.881 

0.980 

0.775 

1.-853 
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where  Pf  is  the  pressure  at  the  detonation  front.  Pf  is  obtained 
by  using  P^  =10poUf  D.  Uf  is  the  particle  velocity  at  the  detonation 
front  and  xs  obtained  from  the  unreacted  Hugoniot.  For  TNT  at  Pq 
= 1.614  g/cc,  the  unreacted  Hugoniot®  is  U = 2.39  + 2.05  u where  U 
is  the  shock  velocity  in  mm/ysec.  The  third  equation  is 


E 


Po"^ 


(19) 


The  integral  in  Equation  (19)  is  evaluated  by  assuming  that  P varies 
linearly  with  t from  P^  to  Pqj.  This  assumption  was  used  because 
the  measured  u(t)  profiles  were  essentially  linear.®  E is  the  energy 
flux  along  the  axis  of  the  charge. 


The  results  of  Equation  (17)  , (18) , and  (19)  are  given  in  Table 
3.  Except  for  data  close  to  the  failure  diameter  the  individual 
results  for  Equations  (18)  and  (19)  are  within  ~10%  of  the  average 
for  each  equation.  This  is  exceptional  considering  the  possible 
error  in  the  measurements.  The  results  for  Equation  (17)  are  not 
quite  as  good.  The  individual  results  are  within  -30%  of  the  average. 
This  is  not  bad  considering  that  there  is  more  error  in  measured  values 
of  P^j  than  D (and  hence  Pf) . Again,  if  an  equation  of  state  existed 
for  the  reaction  products  and  intermediate  products,  then  P^;j  could  be 
calculated  from  and  perhaps  be  more  accurate  than  the  measured  values. 
Thus,  the  indications  are  that  t can  reasonably  be  related  to  P^j  as 
follows : 

T = A"  V(Pcj/PcJo)^  (20) 


where  Equation  (20)  has  been  normalized  to  PqJo* 


An  equation  involving  only  ^CJ  and  R can  now  be  obtained. 
Substituting  Equation  (3)  into  Equation  (11)  results  in 


(21) 


Substituting  Equation  (15)  for  c and  Equation  (20)  for  x there  results 


6.  Dremin,  Saviov,  Trofimov,  and  Shvedov,  Detonation  Waves  in 
Condensed  Media,  (English  Translation) . 

8.  Coleburn  and  Liddiard,  J . Chem.  Phys . , 44 , 1929  (1966). 
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R 

mm 


Equation  (17) 


Equation  (18)  Equation  (19) 


2 2 , 

Kbar  -cc-ysec  /mm-gm 


= 1.62  g/ cc 


30 

7.7  X 10^ 

1.43  X 10^ 

1.065  X 10 

20 

8.0 

1.55 

1.123 

14 

7.7 

1.61 

1.130 

12.5 

6.4 

1.61 

1.058 

11 

5.6 

1.49 

0.959 

10 

5.6 

1.57 

0.997 

9 

5.0 

1.42 

0.895 

3 

4.5 

1.40 

0.854 

Po  = 

1.60  g/cc 

30 

2 

8.4  X 10 

1.55  X 10^ 

1.139  X 10 

20 

7.4 

1.51 

1.082 

17.5 

6.7 

1.45 

1.015 

15 

6.6 

1.31 

0.937 

13.75 

5.0 

1.19 

0.800 
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1/2 

fe)  ■ 


I * II » t I I 


= 1 - 


A'A"A 


fe) 


3/2 


or 


CJ 

XJo 


) 


1/2 


= 1 - 


3/2 


(22) 


(23) 


where  B = A'A"A’ ' ' . 

FAILURE  DIAMETER  PREDICTION 

The  main  question  is:  Does  Equation  (23)  predict  a point  of 

failure?  And  if  so,  how  does  it  compare  with  experiment?  The  first 
question  is  the  same  as  asking;  Does  Equation  (23)  have  an  infinite 
derivative  in  the  CPcj/^CJo) ^ versus  1/R  plane?  To  determine  this 
the  derivative  d (Pcj/PCJo) ^/^/d (1/R)  must  be  obtained.  For  simplicity 
let  (Pcj/Prio) ^ ^ and  1/R  = q.  Then  Equation  (23)  becomes 


= 1 - ^ 
x-’ 


(24) 


or 


x*^  - x^  = - Bq 


:25) 


Differentiating  Equation  (25)  yields 


4x^dx  - 3x^dx  = - Bdq 


(26) 


or 


dx 

dq 


- B 

4x^-3x^ 


(27) 


dx  ^ - B/4 

dq 


At  failure 


x‘(x-3/4) 


(28) 


dx 

dq 


“ x^  = 0 . 75  • 
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Thus  Equation  (23)  has  a vertical  node  in  the  versus  1/R  plane 

when 

1/2 

= 0.75  (29) 

failure 


= 0.56  . (30) 

failure 


Thus  Equation  (23)  predicts  that  there  is  value  of  R,  designated  R^, 
for  which  no  solutions  exist  when  R<Rf.  Therefore,  Rf  is  the  predicted 
failure  radius  of  an  explosive. 

The  answer  to  the  question  on  comparison  with  experimental  results 
can  be  obtained  for  cast  TNT,  at  least.  First  of  all.  Equation  (23) 
predicts  that  the  radius  will  be  reached  when  Prj/Prjo  “ 0.56.  As 
was  mentioned  earlier,  failure  occurs  at  Pcj/^CJo  “ 0.56  and  0.60 

for  1.62  g/cc  and  1.60  g/cc  cast  tnT,  respectively.  Thus  as  far  as  the 
ratio  of  ^cj/^CJo  concerned.  Equation  (23)  works  well  for  cast  TNT. 

1/2 

A calculated  curve  in  the  (P^j/Pp^  ) versus  1/R  plane  can  be 
obtained  as  follows.  A value  for  6 is  obtained  by  using  one 

(Pgj/Pcjp)  ^/^ , ^ datiim  point  and  substituting  it  into  Equation  (23)  . 
Using  this  6/  the  curve  is  generated  by  varying  either  (Pcj/^CJo^ 
or  R and  solving  for  the  other  parameter.  Figure  1 presents  the 
calculated  and  experimental  results  for  1.62  g/cc  TNT.  The  solid 
circles  are  from  Table  2 and  were  calculated  using  Dremin's  data.® 

The  solid  curve  was  calculated  using  Equation  (23)  and  the  datum  point 
at  R = 30  mm  with  6 = 1.176  mm.  The  predicted  Rf  associated  with  this 
curve  is  11.5  mm  while  the  measured  failure  radius  is  8 mm.  The 
predicted  R^  is  44%  too  large.  Note,  however,  that  the  solid  curve  in 
Figure  1 passes  through  only  the  daltum  point  used  to  generate  it.  All 
the  other  data  points  lie  to  the  right  of  the  solid  curve.  The  dashed 
curve  in  Figure  1 was  generated  using  the  datum  point  at  R = 20  mm 
(S  = 0.900  mm).  This  curve  appears  to  fit  the  data  better  and  Rf  is 
8.5  mm  or  6.3%  larger  than  the  measured  value.  The  difference  in  the 
calculated  curves  is  due  to  the  error  in  the  CJ  pressure,  either  the 
measured  Pcj,  or  the  calculated  Pcjo'  both.  Dremin^  unfortunately 
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Figure  2 presents  the  calculated  and  experimental  results  for 
1.60  g/cc  TNT,  the  solid  circles  are  from  Table  2.  The  solid  curve 
was  generated  using  Equation  (23)  and  the  datum  point  at  R = 30  mm 
with  6 = 1.475  mm.  The  predicted  R,  is  13.98  mm  compared  to  the 
measured  value  of  13.75  or  1.7%  too^large.  This  curve  fits  the  data 
so  well  that  no  other  points  were  used  for  comparison. 

In  Table  2,  the  variation  of  A'  with  (Pcj/^CJo^ presented. 

A'  varied  by  a factor  of  1.46  for  1.62  g/cc  and  1.73  for  1.60  g/cc 
TNT.  In  Table  4,  the  calculated  values  of  6 for  each  datum  points  are 
presented.  For  1.62  g/cc  TNT,  8 varies  by  a factor  of  1.39;  for  1.60 
g/cc  TNT,  8 varies  by  1.08.  For  1.62  g/cc  TNT,  the  variation  of  8 
and  A’  are  essentially  the  same.  For  1.60  g/cc  TNT,  8 is  essentially 
constant  while  A'  varies  by  73%.  Note  also  that  the  actual  values  of 
S and  A'  are  very  close.  However,  until  the  physical  meaning  of  8 
is  explained,  no  physical  significance  is  attached  to  this. 

PHYSICAL  INTERPRETATION  OF  THE  RESULTS 

A theory  is  presented  which  relates  the  Chapman- Jouquet  pressure 
of  a detonating  explosive  to  its  charge  radius  and  also  predicts  the 
failure  radius.  The  question  remains;  What  is  the  physical  meaning 
of  this  theory?  First,  this  theory  was  developed  by  comparing  the  ratio 
of  the  total  energy  available  to  the  actual  energy  which  is  unaffected 
by  lateral  rarefaction.  This  resulted  in  Equation  (11) . The  same 
equation  results  if  the  ratio  of  the  energy  flux  through  the  unaffected 
C-J  surface  to  the  energy  fJux  through  the  detonation  front  is  used 
as  a starting  points.  The  important  point  about  either  approach  are; 

(1)  the  unaffected  energy  is  the  governing  parameter,  and  (2)  the 
Chapman-Jouquet  pressure  is  the  measurable  parameter  most  sensitive  to 
changes  in  energy. 

The  second  crucial  step  was  the  assumption  that  P^t  = constant 
holds  in  the  reaction  zone  over  the  range  of  pressure  considered.  An 
alternate  way  of  stating  this  is  that  the  reaction  along  a streamline, 
which  is  not  affected  by  lateral  rarefaction,  will  proceed  until  all 
the  energy  is  released.  The  time  required  for  this  reaction  process 
is,  in  turn,  directly  dependent  on  the  pressure. 

SUMMARY  AND  FUTURE  WORK 

A one  parameter  equation  (Equation  (23))  has  been  developed 
which  relates  Pcjo/  This  equation  fits  the  experimental 

data  for  cast  TNT  reasonably  well  when  the  parameter  6 is  obtained 
using  one  datum  point  and  Pcjq-  This  equation  also  predicts  the 
failure  diameter  of  cast  TNT  reasonably  well  when  the  error  in  the 
pressure  is  taken  into  account.  Unfortunately,  no  complete  set  of 
data  containing  Pqj,  R,  and  the  failure  diameter  has  been  found  for 
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Table  4 

Values  of  6 for  Cast  TNT 


= 1.62  g/cc 


30 

0.955 

1.176 

20 

0.947 

0.900 

14 

0.919 

0.880 

12.5 

0.866 

1.088 

11 

0.845 

1.029 

10 

0.815 

1.001 

9 

0.795 

0.927 

8 

0.751 

0.844 

Average 

=0.981 

= 1.60 

g/cc 

30 

0.941 

1.475 

20 

0.897 

1.487 

17.5 

0.887 

1.380 

15 

0.842 

1.415 

13.75 

0.775 

1.440 

Average  = 1.439 
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a non-ideal  explosive.  Future  experimental  work  will  obtain  such 
a set  of  data  to  check  Equation  (23)  for  a non-ideal  explosive. 

The  derivation  of  Equation  (23)  relies  heavily  on  the  assumption 
that  P^j  T is  a constant,  Equation  (20) . The  data  for  cast  TNT 
indicates  that  this  assumption  is  reasonable.  Future  theoretical  work 
will  investigate  this  relationship. 
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FIG.  1 COMPARISON  OF  PREDICTED  AND  EXPERIMENTAL  P^  VS.  1/R  RESULTS  FOR  CAST  TNT,  Po  = 1.62  g/cc 
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FIG.  2 COMPARISON  OF  PREDICTED  AND  EXPERIMENTAL  P'^  VS.  1/R  RESULTS  FOR  CAST  TNT,  = 1.60  9/cc 


NSWC/WOL/TR  76-158 


DISTRIBUTION  LIST 


Copies 


Chief  of  Naval  Material 
Washington,  D.C.  20360 
MAT-03T 


Commander 
Naval  Air  Systems 
Department  of  the 
Washington,  D.C. 
AIR-350 
AIR-330 


Command 

Navy 

20361 


Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  D.C.  20362 
SEA-0332 
SEA-0  3 3 3A 
SEA-09G32 
Director 

Strategic  Systems  Project  Office 
Department  of  the  Navy 
Washington,  D.C.  20376 
SP-273,  R.  M.  Kinert 
SP-27311,  E.  L.  Throckmorton,  Jr. 


Commander 

Naval  Weapons  Center 
China  Lake,  California  93555 
D.  Lind 

Technical  Library 
H.  D.  Mallory 
G.  Green 


Director 

Naval  Research  Laboratory 
Washington,  D.C.  20375 

Technical  Information  Section 


Director 

Office  of  the  Secretary  of  Defense 
Advanced  Research  Projects,  Agency 
Washington,  D.C.  20301 

Commanding  Officer 
Naval  Weapons  Station 
Yorktown,  Virginia  23691 
R&D  Division 
Code  50 


1 


2 

2 


2 


i 


NSWC/WOL/TR  76-158 


DISTRIBUTION  LIST  (Continued) 


Commanding  Officer 
Naval  Propellant  Plant 
Indian  Head,  Maryland  20640 
Technical  Library 

Commanding  Officer 
Naval  Weapons  Evaluation  Facility 
Kirtland  Air  Force  Base 
Albuquerque,  New  Mexico  87117 
Code  AT-7 

Superintendent  Naval  Academy 
Annapolis,  Maryland  21402 
Library 

Hercules  Incorporated 
Allegany  Ballistics  Laboratory 
P.O.  Box  210 

Cumberland,  Maryland  21502 
Library 

AMCRD 

5001  Eisenhower  Avenue 
Alexandria,  Virginia  22302 

Redstone  Scientific  Information  Center 
U.S.  Army  Missile  Command 
Redstone  Arsenal,  Alabama  35809 
Chief,  Documents 

Commanding  Officer 
Picatinny  Arsenal 
Dover,  New  Jersey  07801 
SARPA-TS-S  #59 
J.  Hershkowitz 

Commanding  General 
Aberdeen  Proving  Ground 
Maryland  21005  • 

BRL 

Commanding  Officer 
Harry  Diamond  Laboratory 
2800  Powder  Mill  Road 
Adelphi,  Maryland  20783 
Library 


2 


Copies 


1 


NSWC/WOL/TR  76-158 


DISTRIBUTION  LIST  (Continued) 


Copies 


Armament  Development  and  Test  Center 
DLOSL/Technical  Library 
Eglin  Air  Force  Base 
Florida  32542 

Energy  Research  and  Development  Administration 
Washington,  D.C.  20545 
DMA 

Director  . • 

Defense  Nuclear  Agency 
Washington,  D.C.  20305 

Research  Director 
Pittsburgh  Mining  and  Safety 
Research  Center 
Bureau  of  Mines 
4300  Forbes  Avenue 
Pittsburgh,  Pennsylvania  15213 
R.  W.  Van  Dolah 

Lawrence  Livermore  Laboratory 
University  of  California 
P.O.  Box  803 

Livermore,  California  94551 
M.  Finger 
E.  James 
£.  Lee 
J . Kury 

Sandia  Laboratories 
P.O.  Box  5800 

Albuquerque,  New  Mexico  87115 
D.  Hayes 

Director 

Los  Aleunos  Scientific  Laboratory 
P.O.  Box  1663 

Los  Alamos,  New  Mexico  87544 
Library 
L.  C.  Smith 
B.  G.  Craig 
A.  Popolato 


DDESB 


Forrestal  Building 
Room  GS  270 
Washington,  D.C. 


20314 


3 


NSWC/WOL/TR  76-158 


DISTRIBUTION  LIST  (Continued) 


Shock  Hydrodynamics,  Inc. 

15010  Ventura  Boulevard 
Snerman  Oaks,  California  91403 
Lr.  L.  Zernow 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  California  94025 
D.  Curran 
M.  Cowpertnwaite 

R.  Stresau  Laboratory,  Inc. 

Star  Route 

Spooner,  Wisconsin  54801 
Rohm  and  Haas' 

Huntsville  Defense  Contract  Office 
723-A  Arcadia  Circle 
Huntsville,  Alabama  35801 
K.  M.  Schuey 

Princeton  University 

Department  of  Aerospace  and  Mechanical 
Sciences 

Princeton,  New  Jersey  08540 
M.  Summerfield 

Chief 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 
Code  473,  R.  Miller 
Code  473,  J.  H.  Satkowski 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  2231^ 


4 


Copies 


